arXiv: 1505.08152vl [cond-mat.quant-gas] 29 May 2015 


Probing the Excitations of a Lieb-Liniger Gas from Weak to Strong Coupling 


F. Meinert,' M. Panfil,^ M. J. Mark,'’^ K. Lauber,' J.-S. Caux/ and H.-C. Nagerl' 

^ Institut fiir Experimentalphysik und Zentrum fur Quantenphysik, Universitdt Innsbruck, 6020 Innsbruck, Austria 
^SISSA-International School for Advanced Studies and INFN, Sezione di Trieste, 34136 Trieste, Italy 
^Institut fur Quantenoptik und Quanteninformation, Osterreichische Akademie der Wissenschaften, 6020 Innsbruck, Austria 
Institute for Theoretical Physics, University of Amsterdam, 1090 GL Amsterdam, The Netherlands 

(Dated: June 1, 2015) 


We probe the excitation spectrum of an ultracold one-dimensional Bose gas of Cesium atoms with repulsive 
contact interaction that we tune from the weakly to the strongly interacting regime via a magnetic Feshbach 
resonance. The dynamical structure factor, experimentally obtained using Bragg spectroscopy, is compared to 
integrability-based calculations valid at arbitrary interactions and finite temperatures. Our results unequivocally 
underly the fact that hole-like excitations, which have no counterpart in higher dimensions, actively shape the 
dynamical response of the gas. 

PACS numbers: 37.10.Jk, 03.75.Dg, 67.85.Hj, 03.75.Gg 


Interacting quantum systems confined to a one-dimensional 
(ID) geometry display qualitatively different behavior com¬ 
pared to their higher-dimensional counterparts H]. Systems 
of strongly interacting electrons recently realized in electronic 
nanostructure devices Elia have evidenced the breakdown of 
Landau’s Fermi liquid theory of quasi-particles in ID, a world 
in which new types of excitations emerge out of the inevitably 
collective nature of the dynamics. The understanding of these 
requires approaches going beyond Landau’s paradigm. The 
best-known, valid for sufficiently small temperatures and en¬ 
ergies, is the Luttinger Liquid (LL) formalism 0. When 
probing dynamical correlation functions, one however typi¬ 
cally leaves this low-energy and large-wavelength limit and 
enters a regime where even recent extensions of the LL for¬ 
malism to higher energies 13161 cannot capture all features. 
Instead, one must rely on nonperturbative calculations to un¬ 
derstand the correct basis of excitations and quantitatively ex¬ 
plain experiments, a recent example being spinon dynamics 
in quantum spin chains da. These systems however lack 
the tunability required to track the whole transformation oc¬ 
curring between the limits of weak and strong coupling. 

Very recently, systems of ultracold bosons have opened up 
new routes to study strong correlation effects in ID d, such 
as the fermionized Tonks-Girardeau (TG) gas of bosons Gol- 
El, primarily due to unprecedented control over system pa¬ 
rameters, e.g. confinement, particle interactions or quantum 
statistics Ha. Moreover, the ID Bose gas with contact in¬ 
teractions is one of the few integrable many-body problems 
that allows for combining experiment with numerically ex¬ 
act studies of the excitation spectrum, making it an important 
cornerstone on the path towards understanding interaction ef¬ 
fects on dynamical correlation functions M- In their semi¬ 
nal work isKia, Lieb and Liniger have shown that next to 
a particle-like mode (Lieb-I mode), which resembles Bogoli- 
ubov excitations in the limit of weak interactions, a second 
mode naturally emerges (Lieb-II mode) that stems from hole¬ 
like excitations in the effective Fermi sea in ID. The coex¬ 
istence of these two types of elementary excitations leads to 
a significant broadening of the dynamical response functions, 



FIG. 1: (color online). Sketch of the experimental setup, (a) A pair 
of retro-reflected laser beams creates an ensemble of ~ 4000 inde¬ 
pendent one-dimensional Bose gases, (b) The excitation spectrum is 
probed by illuminating the gas with a pair of Bragg laser beams. 
(c),(d) Zero-temperature dynamical structure factor S{k, co) (value 
shown in gray scale) for a moderately (c) (y = 3.3) and strongly (d) 
(y = 45) interacting homogeneous gas. The solid (dashed) line in (d) 
shows the dispersion of the Lieb-I (Lieb-II) mode. Insets indicate 
averaging over an ensemble of trapped systems. The thin lines show 
fixed momentum cuts at representative densities. The correspond¬ 
ing values for k/kf are indicated as vertical lines in the k-co-plane. 
The thick line shows the averaged response 5(co) in a local density 
approximation. 


clearly visible in the strongly interacting regime lEllIsl (see 

Fig-0- 

In this Letter, we measure the dynamical structure factor 
(DSF) of the Lieb-Liniger Bose gas realized with ultracold 
atoms confined to ID quantum tubes with widely tunable in¬ 
teractions. The analysis is based on careful disentangling of 
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the experimental traits and allows us to identify the role of the 
Lieb-Liniger dynamics in shaping the response of the system. 
Comparison of the measured spectra with state-of-the-art nu¬ 
merical calculations fnlllli ranging from the weakly to the 
strongly interacting regime allows for a clear distinction be¬ 
tween interaction and temperature effects and demonstrate the 
contribution of the Lieb-11 type excitations to the response of 
the system. 

Our experiment starts with a Cesium Bose-Einstein conden¬ 
sate (BEC) of typically 1.1 x 10^ atoms confined in a crossed 
dipole trap ll20l ISTl . The BEC is adiabatically loaded into 
an array of « 4000 quantum wires created via two mutu¬ 
ally perpendicular retro-reflected laser beams at a wavelength 
X = 1064.5 nm. At the end of the ramp the lattice depth along 
the horizontal direction is Vxj = SOE/;, creating an ensem¬ 
ble of independent one-dimensional ’’tubes” with a transversal 
trap frequency Ox = 27r x 14.5 kHz oriented along the vertical 
z-direction (Eig. [^a)). Here, Er — h^/{2mX^) is the photon 
recoil energy with the mass m of the Cs atom. During lat¬ 
tice loading the scattering length is set to = 173(5)ao 
via a broad Eeshbach resonance ll22l . In the deep lattice we 
then ramp Us within 50 ms to the desired value in the range 
lOflo ^ ^ 900ao to prepare the tubes close to the adiabatic 

ground state. The ramp of Us is carefully adapted to avoid any 
excitation of breathing modes. 

The gas in each tube is described by the Lieb-Liniger 
Hamiltonian Ida 

^ = £ 5{zi-Zj), (1) 

with giD = 2h(0±as {I — 1.0326as /the coupling 
strength in ID ifTTl l23l l24l and a± = sjhj(md)\/) the trans¬ 
verse harmonic oscillator length. The system is conveniently 
described in terms of the dimensionless interaction parameter 
7 = rngi-o/{fP'nxY)), where mid denotes the one-dimensional 
line density 0 . The density sets the characteristic Eermi 
wave-vector kp = Ttnio of the system. In our experimental 
setup we have to consider two sources of inhomogeneity. 
Eirst, the tubes are harmonically confined along the longi¬ 
tudinal direction with a trap frequency (O^ = 27t x 15.8(0.1) 
Hz. This gives rise to an inhomogeneous density distribution 
in each quantum wire. Second, the loading procedure leads 
to a distribution of the number of atoms across the ensemble 
of ID systems fT2\ . Eor comparing measurements with 
theoretical predictions both effects can be accounted for by 
averaging over homogeneous subsystems in a local density 
approximation (EDA) (see insets to Eig. [TJc) and (d)). 

We probe the spectrum of elementary excitations via two- 
photon Bragg spectroscopy ll25l . In brief, the sample is illu¬ 
minated for 5 ms with a pair of phase coherent laser beams at a 
wavelength Ab « 852 nm and detuned by « 200 GHz from the 
Cs D2 line. The beams intersect at an angle (j) at the position 
of the atoms and are aligned such that the wave-vector dif¬ 
ference points along the direction of the tubes (Eig.[TJb)). Its 
magnitude k = 4;r/AB sin( 0 / 2 ) sets the momentum transfer. 
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FIG. 2: (color online). Bragg-excitation spectra for different val¬ 
ues of the ID interaction strength 7 . (a)-(e) Transferred energy 

~ (p^) (normalized to unit area) as a function of the Bragg detun¬ 
ing CO. The scattering length is set to ISaq (a), 173ao (b), 399ao 
(c), 592ao (d), and 819flo (e)> giving an average 7 of 

0.12(5) (9(1)), 3.3(1) (4.4(2)), 11.0(4) (3.7(1)), 21.7(7) (3.5(1)), 
and 45(1) (3.3(1)), respectively. Solid lines show fits to the data us¬ 
ing a Gaussian multiplied by to. The vertical dashed lines indicate 
the position of the Lieb-I (LI) and Lieb-II (L2) mode calculated with 
the averaged values for 7 . The dashed line in (e) shows the calcu¬ 
lated response for ensemble-averaged trapped TG gases, (f) Cen¬ 
tral excitation energy co^ extracted from the Gaussian fit model as a 
function of as (circles). Solid (dashed) lines denote the calculated 
position of the Lieb-I/II modes (Bogoliubov mode). The dotted lines 
indicate the energy of particle and hole excitations in the TG limit, 
(g) £Bfc extracted from the data shown in (a) (triangle), (b) (square), 
(c) (diamond), and (d) (inverted triangle) plotted in the dimension¬ 
less energy-momentum plane. Solid (dashed) lines denote the Lieb-I 
(Lieb-II) mode. The shaded area shows the continuum of excitations 
in the TG limit. Vertical lines in (f) and (g) give the fitted FWHM. 


while a small frequency detuning Cd between the laser beams 
defines the energy transfer to the system. In linear response, 
the energy absorbed from the Bragg lasers for a fixed pulse 
area AE{k,(o) directly relates to the dynamical structure fac¬ 
tor 5(k,a)) = JdxJ c/fe'®'^'^^(p(x,f)p(0,0)) at finite temper¬ 
ature T via A£(k, o) Tico{l — co) with Boltz- 

man’s constant kp ll26ll . 

In our experiment, we probe the ensemble of ID tubes at a 
fixed k — 3.24(3) pm^*, which is comparable to typical mean 
values for kp averaged over the sample ll27ll . The absorbed 
energy as a function of (0 is measured in momentum space. 
For this, we switch off the lattice potential quickly (within 
300ps) and allow for 50 ms time-of-flight at a small positive 
scattering length of « 15ao. From the integrated line density 
along the z-direction of the tubes we extract (p^) and plot it 
as a function of co. The result for five different values of 7 is 
depicted in Fig.|^a)-(e). 

The datasets cover the regime from weak to strong interac¬ 
tions 0.1 < 7 < 50 probed at 0.3 < k/kp < 1. The variation 
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FIG. 3: (color online). Comparison of the Bragg-excitation spectra with theoretical predictions at finite temperature. Symbols depict the 
data presented in Fig.|^for y =3.3(1) (a), 11.0(4) (b), 21.7(7) (c), and 45(1) (d). The vertical dashed lines indicate the position of the Lieb-I 
and Lieb-II modes calculated with the averaged values for y. The solid lines show the ensemble averaged response derived from the finite 
temperature dynamical structure factor, (e)-(h) Reduced analysis of the corresponding experimental data in the top row with theoretical 
predictions for different temperatures. 


in k/kp ensues from the change of the density distribution in 
the tubes with increasing a*, evolving from a Thomas-Fermi 
profile at weak interactions towards the TG profile at strong 
interactions. The values for y and kp given in Fig. [^denote 
the average over the ensemble of ID systems using the mean 
niD in each wire calculated for T = 0 from the solution of the 
Lieb-Liniger integral equations in LDA ll22l . Error bars re¬ 
flect mainly a ±10% uncertainty in the total atom number. A 
clear interaction-induced broadening and shift of the spectra 
with increasing y is observed in accordance with the calcu¬ 
lated position of the Lieb-I and Lieb-II modes (vertical dashed 
lines) ca. We compare the dataset in the strongly interacting 
regime (Lig. [^e)) to the calculated response for an ensemble 
of trapped TG gases at zero temperature averaged over the en¬ 
semble of tubes (dashed line) ||22l|^. The agreement with 
the data underlines the contribution of Lieb’s hole-like excita¬ 
tion to the dynamical response. 

Prior to a detailed discussion on the exact lineshape for fi¬ 
nite y and finite T, we attempt a simplified zero-temperature 
analysis of our spectroscopy signal. A function a)(^(a)c) 
is fit to the data (solid lines), where the DSL averaged over 
the ensemble of tubes is approximated by a simple Gaussian 
function centered at (Oc- The extracted (Oc as a function of 
as is shown in Lig.|^f) (circles). The vertical lines denote the 
fitted full width at half maximum (LWHM). We find the spec¬ 
tral weight of the data within the Lieb-I and Lieb-II modes 
(solid lines) calculated using the ensemble averaged values 
for y. Bogoliubov’s quasi particle energy alone (dashed line) 
does not explain the observed (Oc with increasing y > 3. We 
summarize the measured spectral position and width for four 
selected values of y in the dimensionless energy-momentum 
plane in Lig.|^g), with (Op = Sp/h — hkp/{2m). Comparing 
to the dispersion relation of Lieb’s particle and hole-like ex¬ 
citation, we observe a clear signature for the contribution of 
both branches with increasing interactions, finally approach¬ 
ing the limit of the excitation spectrum of the fully fermion- 
ized Bose gas (shaded area). 


We now turn to a detailed analysis of the spectral lineshape 
as a function of interaction strength and temperature. The 
effect of temperature on the measurement of the DSL arises 
from two distinct contributions. Lirst, the DSL itself is tem¬ 
perature dependent. This becomes most evident in the TG 
limit of fermionized bosons. Lor kpT <^ep the effective Lermi 
sea in quasi-momentum space has a sharp edge giving rise to a 
homogeneous continuum of excitations lying between Lieb’s 
hole- and particle-like modes. When kpT ~ ep the Lermi edge 
washes out, resulting in a smoothening of the DSL with its 
spectral weight being shifted to higher energies am. Second, 
finite temperature affects the density distribution in our ID 
systems leading to a decreasing mean nio with increasing T 
for fixed Ug. This changes the average k/kp at which the tubes 
are probed. 

Lor our theoretical analysis we take both effects into ac¬ 
count. Lirst, we calculate the density distribution in each of 
the tubes at a temperature T by numerically solving the Yang- 
Yang integral equations for the ID Bose gas ll2^ . The DSL 
is evaluated at finite T via the ABACUS algorithm, a Bethe 
Ansatz-based method to compute correlation functions of in- 
tegrable models 1^ . The effect of the trapping potential is 
incorporated by making a LDA for each tube. The response 
of the array of tubes is finally calculated by weighting the con¬ 
tribution of each subsystem by the number of atoms ll22ll . 

The result of our theoretical analysis for four different val¬ 
ues of y is presented in Lig. and compared to the corre¬ 
sponding experimental data (taken from Lig. (b-e)). Al¬ 
though finite temperature leads to small shifts and broadening 
of the excitation spectrum, the most relevant contribution to 
the spectral shape stems from the broadening of the dynam¬ 
ical structure factor with increasing interactions. The analy¬ 
sis underlines the contribution of hole-like excitations to the 
overall response when entering deep into the strongly corre¬ 
lated regime. Lurther, a reduced analysis of our data with 
the computed spectra serves as a thermometry tool in the tubes 
and points to gas temperatures in the range of 5 to 10 nK. A 
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FIG. 4: (color online). Momentum distribution n{p) for resonant 
and off-resonant Bragg excitation measured at weak and strong in¬ 
teractions. (a)-(c) Integrated line density after 50 ms time-of-flight 
when exciting the ID systems at 7 = 0.12(5) (a), 3.3(1) (h), 45(1) 
(c) for red-detuned (left), resonant (middle) and blue-detuned (right) 
excitation with respect to the Bragg resonance. Datapoints show 
the average of four measurements. The frequencies given denote 
the Bragg detuning co/{2n). The red line shows a fit to the central 
part of the line density around p = 0 using a linear combination of 
a Lorentzian and Gaussian profile to extract the width w. The red 
shaded area denotes the fraction of atoms asymmetrically scattered 
to higher momentum states, (d)-(f) Imparted mean momentum (p) 
extracted from n(p) (squares) and relative change in the fitted width 
Sw of the central peak around p = 0 (triangles) as a function of the 
detuning co/{27t) for 7 = 0.12(5) (d), 3.3(1) (e), and 45(1) (f). Solid 
lines are interpolating functions to guide the eye. 


moderate increase in temperature is seen for increasing values 

of7Ea. 

So far, we have characterized the excitations in the gas by 
measuring (p^). Now, we analyze the full momentum dis¬ 
tribution n{p) of the excited ID Bose gases. In Fig. |^a)-(c) 
we plot the atomic line density after time-of-flight, integrated 
transversally to the direction of the tubes, which reflects the 
in-trap momentum distribution of the atomic ensemble. The 
measurements are taken at three different values of 7 , rang¬ 
ing from the weakly to the strongly interacting regime, and at 
Bragg excitation frequencies co slightly below (left column), 
just at (central column), and slightly above (right column) the 
peak of the resonance. 

First, we recognize a dramatic qualitative change in n{p) 
with increasing 7 for all detunings presented. In the weakly 
interacting regime (Fig. Qa)), we observe a clear particle¬ 
like excitation located at p — hk as expected from the non¬ 
interacting limit. Yet, with increasing interactions this feature 
smears out and evolves towards an overall broadening of n{p), 
indicative of a strong collective response of the system to the 
Bragg pulse. This observation demonstrates one of the key 


features of ID systems; any excitation to the system is neces¬ 
sarily collective, and therefore leads to energy-broadened re¬ 
sponse functions, in contrast to sharp coherent single-particle 
modes. This broadening however only becomes clearly visi¬ 
ble for strong enough interactions, where the hole-like modes 
become dynamically relevant. 

In a further measurement, we attempt to quantify the re¬ 
sponse in momentum space in more detail. Note that our pre¬ 
vious measurement of {p^) captures both a broadening of the 
momentum distribution as well as an increase in the mean mo¬ 
mentum {p). In order to separate both contributions spectro¬ 
scopically, we plot the relative change of the width w of 
the central part of n{p) around p = 0 and (p) as a function of 
CO for different values of 7 in Fig. |2d)-(f). The data indicate 
how collective excitations in the gas, expressed via energy de¬ 
position in w rather than in (p), become dominant with in¬ 
creasing 7 . Interestingly, the two curves peak at different val¬ 
ues for CO. This observation is confirmed by the momentum- 
space character of elementary excitations in the interacting 
ID systems, which changes from a collective broadening for 
the Lieb-II mode to a particle-like feature for the Lieb-I mode 
with increasing co ll 22 l . 

In conclusion, we have measured the excitation spectrum of 
a strongly correlated ID system for a wide range of interac¬ 
tion parameters. Comparison with integrability-based calcu¬ 
lations at finite temperature allows for a direct observation of 
the contribution of the collective Lieb-II mode to the DSF. Our 
results demonstrate the successful application of an integrable 
model to analyze dynamics of the ID Bose gas. Furthermore, 
the collective nature of elementary excitations in ID with in¬ 
creasing interaction strength has been demonstrated through 
an analysis of the momentum distribution. Our results pose 
questions on the time evolution and potential equilibration of 
these collective excitations. This could be seen as an alterna¬ 
tive quantum cradle setting 1311 in which, instead of colliding 
two highly energetic clouds of atoms, relatively low energy 
excitations propagate through the system, whose individual 
features can then be more easily studied. 
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FIG. 6 : ID line-density distribution of a tube with N = 30 atoms at 
flj = 150ao (a) and flj = 800ao (b)- The blue line shows the result 
from the numerical solution of the Lieh-Liniger system. The dashed 
(dotted) gray lines depict the analytic results for the profile in the 
Thomas-Fermi (Tonks-Girardeau) limit m Here, = 2;r x 15 kHz 
and cOj = 15 Hz. Finite temperature profiles using the Yang-Yang 
thermodynamic equations are shown for T = lOnK (orange line) and 
T = 30nK (red line). 


FIG. 5: (a) Atom number distribution over the ID tubes for a total 
atom number A = 1.1 x 10^ in the BEG and a 3D scattering length 
during loading 030 = 173flo- (b) Number of tubes that are filled with 
Nij = N' atoms, (c) Total number of atoms in tubes filled with Nij = 
N' atoms. 


SUPPLEMENTARY MATERIAL: PROBING THE 

EXCITATIONS OE A LIEB-LINIGER GAS EROM WEAK 
TO STRONG COUPLING 

Lattice depth calibration and error bars 

The lattice depth Vx,y is calibrated via Kapitza-Dirac 
diffraction. The statistical error for Vx^y is 1%, though the sys¬ 
tematic error can reach up to 5%. 

The scattering length is calculated via its dependence on 
the magnetic field m with an estimated uncertainty of ±5ao 
arising from systematics in the magnetic field calibration and 
conversion accuracy. Additionally, the magnetic field gradient 
for sample levitation nzm leads to a variation of less than 
±3ao across the sample. 


that interactions are sufficiently small during the loading pro¬ 
cess so that all tubes are in the ID Thomas-Fermi (TF) regime, 
the local chemical potential in each tube reads 

M/.j = M {0)1 + CO^f) , 

with m the mass of the Cs atom and A = 1064.5 nm the wave¬ 
length of the lattice light. From jj^ij the atom number in each 
tube can be derived via 


^Uj 


UV2 




2/3 


with the ID coupling strength 


giD = 2ncoj_as 


1-1.0326 — 


-1 


Here, a± = y/h/{m(0±) is the radial oscillator length. The 
global chemical potential is calculated iteratively from the 
condition N = 


Atom number distribution across tbe ID tubes 

In order to calculate a mean value of the interaction param¬ 
eter 7 and the Fermi wave vector kp we need to know the atom 
number distribution across the ensemble of tubes, as shown in 
Fig .|^ We follow the calculation presented in lU. 

The EEC is adiabatically loaded from a crossed dipole trap 
into the optical lattice. When the lattice is fully ramped 
up to its final depth Vx^y = 3 QEr, the laser beams forming 
the array of tubes give rise to an additional background har¬ 
monic confinement. The total background harmonic confine¬ 
ment of lattice and dipole trap laser beams is measured to 
(Ox = 2nx 15.4(0.1)Hz, co,, = 2nx 20.1(0.1)Hz, and (o, = 
2k X 15.8(0.1) Hz. We deduce the atom number Nij for the 
tube ((, 7 ) from the global chemical potential pL. Assuming 


Density profile in the tubes 

For r = 0 we model the ID density distribution n(z) in each 
tube individually by numerically solving the Lieb-Liniger sys¬ 
tem and making a local density approximation EUa. For 
r > 0 we solve the Yang-Yang thermodynamic equations of 
the ID Bose gas making a local density approximation to cal¬ 
culate n{z,T) Q. Representative examples for two different 
values of Us at zero and finite temperature are shown in Fig.|^ 

Mean 7 and kp 

From the atom number distribution, we calculate the den¬ 
sity profile for each tube individually as described above. The 
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scattering length (ao) 


FIG. 7: (top row) Relative atom number in tubes with local interac¬ 
tion parameter jij (a) and a local value for k/k‘p (b). As an example, 
the distributions are calculated for = 399ao A = 1.1 x 10^ re¬ 
sulting in a global y{kf) of 11 (3.7^m^*). (bottom row) Global val¬ 
ues for y (c) and kf (d) as a function of the scattering length Aj. The 
points show the calculated values with iiid computed numerically 
within the Lieb-Liniger model using a local density approximation. 
The dashed line in (d) is calculated with the density in the Thomas- 
Fermi limit while the dotted line shows the result using the density 
for a Tonks-Girardeau gas. The solid lines are interpolation functions 
to guide the eye. In (a)-(d) all quantities have been calculated using 
the mean value for aid in each tube. 


mean ID density in each tube nj^ then delivers local and mean 
values for y and kp 


r.j = 


mgiD 

h^n 






and 






F ■ 


1-.I 


As an example, the distribution of atoms over tubes with local 
Yij and k/k'p^ is depicted in Fig. |^a) and (b) for a specific 
value of dj. In Fig.j^c) and (d) we show the mean value of y 
and kp as a function of Us- 


Sampling of the dynamical structure factor over the array of 
tubes and comparison with the experimental data 

In linear response, the energy AE{k,( 0 ,T) dumped into a 
single ID gas after the Bragg pulse relates to the dynamical 
structure factor S{k, ( 0 ,T) via |[8]| 

AE{k, CO, T)ocrico{l- co, T). 


Provided that the system is probed at sufficiently large mo¬ 
mentum k the DSF of the trapped system can be computed in 
a local density approximation 

L 

S{k,co,T) = — J Shom{k,CO,T;n{z))dz, 

-L 

with L the system length and 5hom the DSF calculated for a 
homogeneous system with uniform density. In practice, we 
find that dividing the profile of each tube in « 10 homoge¬ 
neous subsystems approximates the DSF of the trapped gas 
sufficiently well. The DSF in each tube with A,,; atoms obeys 

the/-sum rule IHl, / a)5(A:,a),7’)t/a) . In combination 

with detailed balance S{k, co) = —co) we find 

-foQ 

J coSik,co,T){l- oc Aj. 

0 

This allows us to calculate the dynamical response of the en¬ 
tire ensemble by first normalizing the DSF of each individual 
tube {i,j) by its /-sum, respectively, and then weighting its 
contribution to the total signal by the number of atoms Nij. 
For a direct comparison with the experimental signal, we 
finally normalize the ensemble averaged response to unit area. 

In the experiment, we extract (p^) as a function of co 
from the momentum space distribution obtained after time- 
of-flight. For a direct comparison, each spectrum is normal¬ 
ized to unit area as the overall signal depends on the intensity 
of the Bragg lasers, which we slightly increase for data taken 
at stronger interactions. A changing offset due to an overall 
broadening of the momentum distribution with increasing y is 
subtracted from the data. This offset does not affect the shape 
of the excitation spectrum and stems from the broadening of 
the unperturbed momentum distribution with increasing y. 


The ABACUS algorithm 

The DSF for a homogeneous system is computed using the 
ABACUS algorithm. The computations are performed for a 
finite system of length L with a finite particle number N and 
with periodic boundary conditions. This means that the ex¬ 
perimental situation is recovered only in the thermodynamic 
limit N,L^ °° with a density n=N/L fixed through the local 
density approximation. We have confirmed that the particle 
numbers N chosen for the computations (up to A = 128) are 
large enough for the results to faithfully represent the thermo¬ 
dynamic limit. The only exceptions are the highest temper¬ 
ature curves in Fig. 3 (b)-(d) where some oscillations due to 
the finite size are still visible at high energies. The correlation 
function is evaluated numerically by summing contributions 
according to the Lehmann spectral representation. The spec¬ 
tral sum is infinite and needs to be truncated. The ABACUS 
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Bragg laser power (mW) Bragg laser power (mW) 


FIG. 8: (a)-(c) Bragg excitation as a function of the laser pulse length 
measured in a weakly interacting 3D BEC (a), in ID at Os = ISqq (b)> 
and in ID at flj = 819ao (c). The laser power in both Bragg beams 
is set to O.SmW (circles) and 0.2mW (squares) in (a) and (b), and 
0.5mW (circles) and 0.3mW (squares) in (c). (d),(e) Bragg excitation 
as a function of the laser power in both Bragg beams measured in a 
weakly interacting BEC (circles), in ID at = 15ao (squares), and 
in ID at as = 819ao (diamonds). Here, the laser pulse length is fixed 
to Sms. The solid lines are exponentially damped sinusoids fit to the 
data. The dashed lines are linear fits in (a)-(c) and quadratic fits in (d) 
and (e) to the initial increase denoting the regime of linear response. 


algorithm performs this operation in an efficient way captur¬ 
ing the most relevant contributions. The error caused by the 
truncation is easily tractable by evaluating the /-sum rule and, 
for the presented data, does not exceed 5% of the total spectral 
weight. 


Regime of linear response 

Measuring the DSF via Bragg spectroscopy requires to 
probe the system in the regime of linear response. This we 
tested experimentally by validating that the excitations cre¬ 
ated depend linearly on the pulse length and quadratically on 
the intensity of the Bragg lasers for the range of parameters 
used in the experiment, see Fig. Si). 



FIG. 9: (a) Heating effects during the ramp of the scattering length 
as on the Bragg spectroscopy data. Bragg spectroscopy data taken 
directly after the ramp of the scattering length to 173ao within 50 ms 
(squares) is compared to data taken after ramping to 819ao and back 
within 100 ms (circles), (b) Full width at half maximum (FWHM) 
extracted via the Gaussian fit (see main article) from the Bragg spec¬ 
troscopy data as a function of y. The dashed (solid) arrow denotes an 
estimate for line broadening due to the finite pulse length (quantum 
finite size effect) in the non-interacting limit. 


Heating effects on the excitation spectra 

Besides the detailed theoretical analysis of finite temper¬ 
ature effects on the measured DSF reported in the main ar¬ 
ticle, we have checked in an experiment that heating during 
the ramp of to large values only marginally influences the 
shape of the excitation spectrum when compared to the effect 
of interactions. In Fig. |^a) we show Bragg excitation spec¬ 
tra taken at moderate interaction strength, corresponding to 
the data shown in Fig. 2(b) of the main article. For the two 
datasets shown, we either ramp directly to 173ao (squares) or 
ramp deep into the regime of strong interactions, a^—SlOaoj 
and back (circles) prior to applying the Bragg pulse. 


Interaction-independent spectral width at small y 

In the limit of weak interactions we observe a width of the 
excitation spectra that levels off at a constant value as shown 
in Fig. I^b). We attribute this to the onset of two effects 
that prevent the observation of a 5-like peak as expected in a 
non-trapped homogeneous system for y —0. First, the finite 
length of the Bragg excitation pulse causes Fourier broaden¬ 
ing, which results in an estimated residual width of « 120Hz 
in the non-interacting limit M- Second, the finite size of 
the sample due to the presence of the harmonic trap leads to 
an uncertainty-limited energy width (not accounted for in the 
LDA) that can be estimated to 5co = 2hk/{ma^) for a sam¬ 
ple of non-interacting particles cni. Here, = sJh/{mO}v) 
denotes the quantum length scale of the trap. 


Momentum distribution of Lieb-I and Lieb-II excitations 

In order to illustrate the effects of additional particle-hole 
excitations on the momentum distribution, we consider two 
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FIG. 10: Momentum distribution of two states excited with a single 
(particle-hole) excitation over the ground state for 7 = 64. In the first 
case (a) the momentum is carried hy a hole-like Lieb II excitation, in 
the second (h) hy a particle-like Lieh I excitation. Both states have 
the same total momentum equal iok/kf = i. The calculation is done 
for a system of 32 particles. The shaded areas indicate the difference 
to the ground-state momentum distribution. 


situations in which the total momentum of an excited state 
(obtained by adding a single particle-hole excitation with total 
momentum equal to the Fermi momentum on the ground state) 
is carried either by the hole (Lieb type II) mode (Fig. [T0|^a)) 
or by the particle (Lieb type I) mode (Fig.[T0|b)). The compu¬ 
tations are performed again using the ABACUS method im. 
The results show that while in both cases the momentum dis¬ 
tribution function is augmented by broadened peaks, the hole 
excitation leads to an overall broadening of n{k) while the par¬ 


ticle excitation appears as a more distinguishable peak. This 
highlights the collective nature of the Lieb-II hole-like excita¬ 
tions and qualitatively reproduces the features seen in Fig. 4 
of the main text. 
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